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1. Introduction 

Elastic deformation is one of the most important con- 
siderations in structural applications of solid materials. 
Indeed, elastic properties are commonly required in 
computer aided design and manufacturing techniques 
to simulate a product's behavior under variable condi- 
tions of stress and temperature. Under such conditions, 
it is desirable to have a means of estimating the value 
of a property continuously at any temperature or stress 
according to the local operating conditions. Tabulated 
data sets, however, are discrete and may be relatively 
sparse, particularly with respect to the dependence on 
microstructure. While interpolation techniques can be 
used with tabulated data if sufficiently extensive data 
tables are available, such approaches are relatively 
cumbersome. A more succinct and efficient approach is 
to use semiempirical analytical models that incorporate 
both material and environmental factors within the 
model. 

An opportunity to construct an analytical representa- 
tion of the elastic moduli data evolved recently from an 
extensive compilation (NISTIR 6853) of the elastic 
property data for polycrystalline oxide ceramics [1]. In 
that work, data were collected from the technical liter- 
ature, either as reported in textual or tabular formats or 



as digitized from graphical formats. Special attention 
was given to the dependence of the moduli on both 
porosity and temperature. 

In the present work, we report the construction of 
analytical representations of the elastic moduli data 
using a single model in which the effects of porosity (<^) 
and temperature (7) are treated simultaneously. Results 
for this model, applied to the data in NISTIR 6853 to 
the extent that sufficient data were available to evaluate 
the parameters in the model, are presented for 24 mate- 
rial specifications. 



2. Model 

To construct a suitable model, we proceed heuristi- 
cally, beginning with the assumption that a separation 
of variables may be applied to the dependence of elas- 
tic moduli on temperature and porosity. For any modu- 
lus, M{T,(p), of a given material composition, it is 
assumed that (p and Tmay be taken as independent vari- 
ables, and hence that we may consider 



M(T,<l>)=M,(T)M^(<l>) 



(1) 
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such that our task is to find suitable representations for 
Mj(T) and M^ (({>). 

2.1 Temperature Dependence 

Empirically, the temperature dependence of Young's 
elastic modulus for most ceramics is relatively simple, 
generally decreasing monotonically with increasing 
temperature. At very low temperature, the slope of the 
modulus with respect to temperature must approach 
zero. On the basis of lattice dynamics. Bom and Huang 
[2] estimated that the elastic constants should vary as P 
at low temperature. Above room temperature, the mod- 
uli generally decrease linearly with increasing temper- 
ature. To describe the behavior from low to high tem- 
perature, Wachtman et al. [3] suggested the empirical 
relation 



E^iT) = E,-bTexpi-T,/r) 



(2) 



in which Eq is Young's modulus at absolute zero, and b 
and 7^0 are parameters to be determined numerically 
from the observed data. Anderson [4] later provided a 
justification of an expression of this form for the bulk 
modulus and noted that the elastic modulus would be 
approximately of the same form if the temperature 
dependence of Poisson's ratio could be ignored. 

Empirically graphs of elastic moduli data vs temper- 
ature exhibit very little curvature except at very low 
temperature. This lack of curvature causes numerical 
fitting routines to be rather insensitive to the exponen- 
tial factor in Eq. (2). Consequently the uncertainty in 
the value of the parameter, T^q, is unacceptably large for 
most of the data used in the present work. For the pres- 
ent purpose, therefore, it suffices to consider only the 
simplified linear model 



M,iT) = M,iO)il-a^r) 



(3) 



with the parameters rewritten as Mj(0) and au for each 
modulus M. 

2.2 Porosity Dependence 

The porosity dependence of the elastic properties of 
solids has been the subject of extensive investigation 
for decades. Numerous studies have examined the role 
of pores as the second component of two-phase solid 
media [5-10]. Those works generally involve an analy- 
sis of the strain field in the composite body under the 
application of an external stress. Alternatively several 
studies [11-17] have observed that stress internally is 



transmitted only over the areas of contact between the 
constituent particles or grains. As the body is densified, 
the contact area increases while the porosity decreases. 
Consequently the porosity dependence of the elastic 
moduli should be governed by the contact area. More 
recently detailed analyses of the effects of pore size 
and pore shape have begun to be performed in finite 
element computer simulation calculations [18,19]. 

In addition to these microstructural modeling efforts, 
many semiempirical analytical models have been pro- 
posed [20-29] and applied [30-37] to represent the gen- 
eral trend of elastic moduli with porosity. Analytical 
models are of considerable interest because of their 
potential use as smoothing and interpolation fiinctions. 
Since these models only relate bulk elastic properties to 
the mean porosity, they generally do not represent 
detailed microstructural effects arising from varying 
pore shape, anisotropy or nonuniformity Their impor- 
tance rests in their capacity to provide highly effective 
descriptions of the trends of the mean properties and 
characteristics of porous media. 

Empirically a simple linear model [20] may be ade- 
quate at very small porosity, but for most brittle mate- 
rials, the elastic moduli vary approximately exponen- 
tially [22] for porosity up to about 30 %. At higher 
porosity, the elastic moduli may deviate significantly 
from an exponential dependence [38]. Several models 
treat porous media as a special case of a two-phase 
medium in which the second phase consists of pores 
[36]. Those models often express the moduh of porous 
materials as ratios, P^{<p)IP2{<p), of polynomials {P^ and 
P2) in the volume fraction of porosity (0). Budiansky's 
self-consistent model [30] is of this type and results in 
a pair of coupled equations for the bulk and shear mod- 
uli. Those relations are explicitly linear in porosity and 
implicitly nonlinear through the self-consistent depend- 
ence on Poisson's ratio, v, which is itself dependent on 
porosity. 

At very high porosity, other issues must be consid- 
ered in determining the influence of porosity on elastic 
moduli. It is self-evident that the volume fraction of 
porosity of a solid material must be less than one 
(0< 1) because the condition 0= 1 corresponds to no 
material at all. As the limit = 1 is approached, the 
contiguity of the assemblage of components becomes 
an important issue since the integrity of an elastic medi- 
um is dependent on the transitivity of forces between 
adjacent material components. Indeed, in studies apply- 
ing percolation theory, analyses of minimum solid 
areas of idealized stackings, and other models focused 
on the stacking of geometric shapes, there arises the 
possibility of a critical porosity, 0^? at which the moduli 
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must vanish [11]. Such studies pertain to the very 
important issue of the validity of interpreting such an 
assembly of material components as an elastic continu- 
um. Phani and Niyogi [26] suggested that if we are to 
allow for a vanishing modulus, then Young's modulus, 
E, should be proportional to a power of (1 - <^/<^c)- 

In the present work, elasticity, as a bulk concept, is 
taken to mean a priori that the spatial connectivity is 
sufficient to allow the bulk material to sustain an 
applied stress. For any such material, without excep- 
tion, the elastic modulus does not vanish. 

Assuming material contiguity, Wagh et al. [27] con- 
sidered a model in which the material was assumed to 
be composed of a network of material chains and inter- 
posed with channels of open pores. For a one dimen- 
sional system, they obtained the closed form expression 



E=E,{\-(i)r 



(4) 



where E is Young's modulus, and E^ and n are 
adjustable parameters. They then used numerical solu- 
tions to verify that the same expression should be valid 
also for a three dimensional system. That conclusion 
was consistent with the results of Gibson and Ashby 
[37] who obtained Eq. (4) for the specific case of cellu- 
lar ceramics, with n = 2 for open cell structures and 
n = 3 for closed cells. 

Among these various models, it may be noted that 
the suitability of the various analytical forms is not 
sharply distinguished over the observed range of poros- 
ity for polycrystalline ceramics. No one model seems to 
have a stronger theoretical justification than the others, 
and the empirical fits to the data are not sharply differ- 
ent. Additionally the general trends of the elastic mod- 
uli data vs porosity, for polycrystalline ceramics, do not 
seem to depend greatly on the nature of the porosity 
since results for specimens from multiple sources con- 
form to a single trend line. Neglecting such details, it is 
possible to derive [39] a simple effective medium theo- 
ry for the porosity dependence of bulk moduli. In this 
approach, the classical model of an ionic solid [40] was 
taken as an idealized, pore free, reference system. That 
choice had the particular virtue of providing a closed 
form expression for the bulk modulus. It was noted that 
the introduction of porosity into such a system must 
increase the molar volume of the material, M/p, where 
M is the molecular mass and p is the bulk density. As a 
result, the mean interaction potential at a site must be 
reduced because the mean interparticle distance is 
increased. To account for this relaxation in the model 
system, the length scale was formally renormalized. 
The renormalized system was then related to the porous 



physical system by imposing the consistency condition 
that the equilibrium volume of the renormalized system 
be equal to the sum of the volume at zero porosity and 
the pore volume. The result was the closed form 
expression 



B=B,{\-<i>r 



(5) 



In this model, the exponent, m, was determined by the 
effective attractive component of the interaction poten- 
tial and can be different from the exponent, n, found in 
the similar expression, Eq. (4), for Young's modulus. 

2.3 The General Model 

The elastic properties of polycrystalline ceramics 
usually are approximately isotropic because of the ran- 
domness of the grain orientations, even when the indi- 
vidual grains are anisotropic. An exception to this usual 
circumstance occurs for textured materials in which the 
microstructure has partially aligned grain orientations. 
In the present work, we consider only polycrystalline 
ceramics that may be treated as isotropic materials. For 
this case, the elastic properties are fully described by 
any two of the elastic moduli. 

Upon viewing the dependence on temperature and 
porosity separately, we have seen that the temperature 
dependence may be represented effectively by Eq. (3). 
For the porosity dependence, there are several alterna- 
tives, but only two of the models, Eq. (4) for the elastic 
modulus and Eq. (5) for the bulk modulus, have been 
derived in closed form from theoretical models. 
Combining these models in the manner of Eq. (1), we 
obtain the general model describing the simultaneous 
dependence of E and B on the variables T and <^. 



E{T,(l>)=E,{\-aT){\-(l>y 



(6) 

(7) 



3. Discussion 



The model represented by Eqs. (6) and (7) has been 
applied to the data in NISTIR 6853, and the results are 
given in Table 1 . An illustration of the typical fit of the 
model is given by the results for magnesium aluminate 
spinel [41-46], Fig. 1 and Fig. 2. 

It should be noted that reports of elastic property data 
in the literature most commonly provide results for the 
elastic modulus and the shear modulus, G. The shear 
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Table 1. Parameter values for the fit of the analytical model, Eqs. (6) and (7), for various oxide ceramics. The valid temperature and porosity 
ranges are indicated. The relative expanded uncertainties (coverage factor A: = 2, 95 % confidence limit) for the computed elastic and bulk moduli 
were estimated as 5 %. Brackets, {}, indicate additional approximations were used, as indicated in the footnotes. Mj. = molar mass, pti^e^ = theo- 
retical mass density of the unstressed single crystal at room temperature 



Material 


M, 


Ptheo 


r range 


Porosity 


^0 


a 


n 


^0 


b 


m 




gmol 


gem 


°C 


range 


GPa 


10^ °C"^ 




GPa 


10^ °C"^ 




AI2O3 


101.961 


3.984 


Oto 1000 


to 0.9 


393 


1.33 


3.06 


241 


0.84 


3.33 


Al^Si^Oi/ 


426.052 


3.17 


to 900 


Oto 0.13 


229 


1.17 


3.33 


166 


{1.16} 


3.15 


BeO 


25.012 


3.01 


to 1400 


Oto 0.16 


386 


0.77 


1.96 


350 


1.18 


1.61 


DyzOa^ 


372.998 


8.161 


to 900 


to 0.2 


186 


1.37 


3.81 


144 


{1.37} 


3.52 


Er203 


382.516 


8.651 


Oto 1000 


to 0.2 


179 


1.14 


2.57 


160 


1.14 


3.08 


Gd203 


362.498 


8.348 


to 1400 


to 0.37 


157 


1.46 


2.32 


114 


1.47 


2.19 


Hf02(c,Pr)^'^ 


See^ 


n/a 


Oto 1500 


to 0.09 


251 


1.21 


2.86 


183 


{1.21} 


3.23 


Hf02(c,Tbf^ 
Hf02(c,X)^'^ 


See^ 


n/a 


Oto 1650 


Oto 0.18 


229 


1.41 


1.78 


186 


{1.41} 


2.78 


SeeJ 


n/a 


to 1500 


to 0.38 


{256} 


{1.52} 


{3.01} 


{200} 


{1.70} 


{4.09} 


Hf02(PSH)^'^ 


See^ 


n/a 


Oto 1600 


Oto 0.12 


{263} 


{2.29} 


{3.47} 


{162} 


{2.29} 


{3.45} 


H02O3 


377.859 


8.414 


Oto 1000 


Oto 0.18 


175 


1.08 


2.60 


155 


0.98 


3.43 


LU2O3 


397.932 


9.423 


Oto 1000 


to 0.34 


204 


1.03 


3.12 


161 


0.24 


4.27 


MgAl204 


142.266 


3.572 


Oto 1200 


to 0.38 


278 


1.98 


3.20 


187 


1.97 


3.57 


MgO 


40.304 


3.58 


to 2500 


to 0.26 


310 


1.63 


3.81 


164 


1.23 


2.64 


SC2O3 


137.910 


3.841 


to 1400 


to 0.3 


229 


1.22 


2.97 


148 


0.98 


2.45 


Sm203 


348.718 


7.748 


to 1300 


to 0.38 


150 


2.00 


2.85 


125 


1.73 


3.45 


Th02^ 


264.037 


10.0 


Oto 1200 


to 0.4 


258 


1.68 


3.32 


187 


{1.66} 


4.18 


Ti02'' 


79.866 


4.25 


Oto 1600 


to 0.35 


286 


1.52 


4.99 


{200} 


{2.20} 


{6.57} 


Tm203 


385.867 


8.889 


Oto 1000 


to 0.24 


185 


0.88 


3.07 


147 


1.63 


2.18 


YBa2Cu306 9 


664.594 


6.37 


-268 to 25 


to 0.5 


150 


1.54 


3.70 


69 


1.84 


3.19 


Y2O3 


225.810 


5.03 


to 1600 


to 0.37 


176 


1.37 


2.47 


147 


1.93 


3.27 


Yb203 


394.078 


9.2932 


Oto 1000 


to 0.27 


199 


0.90 


2.61 


155 


1.24 


2.83 


Zr02(m)"^ 


123.223 


5.6 


Oto 1000 


to 0.2 


244 


2.86 


3.79 


170 


3.19 


3.49 


Zr02(c)^'" 


See^ 


n/a 


Oto 1600 


to 0.2 


{227} 


{1.50} 


{2.59} 


{183} 


{1.48} 


{4.31} 



^Neither B{T) nor G{T) was known. Parameters were estimated using Uq = a^. 

^ Parameters estimated using data from specimens with differing dopants. 

^ Optimization routine did not converge. Apparent midrange values were selected manually. 

Monoclinic structure. 
^ Cubic structure, Hf02'xPr203. 
^M, = 210.489 + 329.814X. 
^ Cubic structure, Hf02-xTb203. 
I" Mr = 210.489 + 365.849X. 

' Cubic structure, Hf02-xX203, X = Er, Gd, Pr, Tb, and Y. 
J M, = 210.489 + xM,(X203). 

^ Partially stabilized hafiiia, Hf02-xX203, X = Er, Eu, and Y. 
^ M, = 210.489 + xMr(X203). 
^^ Monoclinic structure. 

" Cubic structure, Zr02'xX203, X = Ca, Pr, Tb, and Y. 
''M,= 123.223 + XMXX2O3). 



modulus for isotropic polycrystalline materials may be 
obtained from E and B as 



3BE 
9B-E 



(8) 



From this relation, it can be seen that G generally will 
not be of the same analytical form as E and B, For 
ceramics, the magnitude of E is typically on the order 



of twice that ofB, Consequently, the relation in Eq. (8) 
can be expanded as 



yielding 



G«-^E 



95 95 



(10) 
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from which it is seen that G may have a different fiinc- 
tional dependence on T and (p, depending on the ratio 
{EI9By 

Similarly, we may note that Poisson's ratio, v, is 
given by 



> 



0.40 
MgAI,0, 

n [41] o [44] 

V [42] O [45] ^0.35 

A [43] O [46] 



\A 



X^----^X 



G 

100 - ^^ 



50 




--^ 



^^." 



A 



0.30 



0.25 



TO 



o 
w 

0.20 -5 

CL 



0.15 
0.10 



0.0 0.1 0.2 0.3 0.4 

Volume Fraction of Porosity (^ 

Fig. 1. Elastic moduli (E, G, and B) and Poisson's ratio (v) of 
MgAl204 vs. porosity, at room temperature. Numbers in square 
brackets, [...], are references for the experimental data. The smooth 
curves comprise the fit of the model, Eqs. (6) and (7). (N.B.: The E 
value from [41] (square point at = 0.23) was treated as an outlier in 
fitting the E data.) 



~6B 



(11) 



and depends directly on the ratio (E/6B). In the present 
work, the magnitudes of the products a^and bTin Eqs. 
(6) and (7) typically were found to have values of about 
0.1 at 1000 °C. Hence, the ratio (E/B) is approximately 



E_ 
B 



B. 



•il-[a-b]T)il-^y- 



(12) 



Consequently Poisson's ratio is not expected to be con- 
stant and may increase or decrease with temperature 
and porosity in a manner that reflects how the depend- 
ence of £' differs from that ofB with respect to the vari- 
ables T and (f). 



4. Conclusion 
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Fig. 2. Elastic moduli of MgAl204 vs. temperature, for various val- 
ues of the porosity. Numbers in square brackets, [...], are references 
for the experimental data. The subscript, "bend," indicates values 
derived from stress-strain bending curves; all other values were 
derived from ultrasonic measurement methods. The smooth curves 
comprise the fit of the model, Eqs. (6) and (7). 



The condensation of a large tabulation of discrete 
data values into a representative analytical model is a 
data evaluation technique that optimizes the utility of 
the collected experiential data. The result is a succinct 
representation that enables the results to be more read- 
ily and consistently integrated into computerized 
design programs and enhances the use of the results in 
distributed data systems. The present work discusses 
the application of that technique to a compilation of 
elastic moduli data for a wide range of polycrystalline 
oxide ceramics. The model used in this work provides 
simultaneous, self-consistent representations of the 
elastic and bulk moduli for polycrystalline ceramics as 
fiinctions of temperature and porosity 
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